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ABSTRACT 


Free electron lasers (FEL's) are devices used worldwide 
for several purposes. In the military, especially in the 
Navy, they can be used for self-defense against missiles, 
and small boats. Installed on a ship, an FEL represents a 
multi-mission, deep magazine, long range weapon. This thesis 
will describe briefly the basic components and principles of 
operation. It also explores, by simulations, the effects of 
changing some of the parameters that generate the laser 
beam. 
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I. 


INTRODUCTION 


A. HISTORY OF THE LASER 

The word "laser" is an acronym for "light amplification 
by the stimulated emission of radiation". Laser is an 
optical source that emits nearly monochromatic light in a 
narrow coherent beam. 

The laser was proposed as a variation of the maser 
("Microwave amplification by the stimulated emission of 
radiation") principle in the late 1950s and was demonstrated 
in July 1960 by Theodore Maiman at Hughes research 
laboratories [1]. 

Albert Einstein, in his 1916 paper, Strahlungs-Emission 
und-absorption nach der Quantentheorie, laid the foundation 
of the laser and the maser. He applied Max Planck's law of 
radiation in terms of the probability of absorption and 
spontaneous and stimulated emission. In 1928, Rudolph W. 
Landenburg confirmed the existence of stimulated emission 
and negative absorption experimentally. In 1939, the 
Russian, Valentin A. Fabrikant, suggested the use of 
stimulated emission to amplify "Short" waves. In 1947, 
Willis E. Lamb and R.C. Retherford demonstrated apparent 
stimulated emission in hydrogen spectra. Finally in 1950 
Alfred Kastler proposed the method of optical pumping, which 
was experimentally confirmed by Brossel, Kastler and Winter 
two years later [1]. 

In 1953, Charles H Townes, James P Gordon and Herbert 
J. Zeiger made the first "Maser" to produce microwave 
radiation. In 1957, Townes and Arthur L. Schawlow began to 
produce emissions in the infrared band and then adapted it 
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to optical wavelengths. In 1958, Prokhorov proposed the 
usage of an open resonator in the maser, and Schawlow and 
Townes settled on an open resonator design[1]. 

The first working laser was made by Theodore Maiman in 
1960 in Malibu, California; a collaboration of several 
institutions. This laser was made with a synthetic ruby 
crystal, pumped by a flashlamp to produce light at 
wavelength 694 nm [1]. 

In 1960, All Javan, William Bennet and Donald Heriot 
made the first gas laser using Helium and Neon [1]. 

The first laser diode was made by Robert N. Hall in 
1962. It was made of gallium arsenide and emitted in the 
near-infrared at 850 nm. The first semiconductor laser with 
visible emission was made the same year by Nick Holonyak, 
Jr. These semiconductor lasers could be used only in pulsed 
operation, and indeed only when cooled to liquid nitrogen 
temperatures (77 K) [1]. 

In 1970, Zhores Alferov in the Soviet Union and Izuo 
Hayashi and Morton Panish of Bell Telephone Laboratories 
independently developed laser diodes continuously operating 
at room temperature, using the heterojunction structure [1]. 

B. HISTORY OF A FREE ELECTRON LASER 

The Free Electron Laser (EEL) appeared in the middle 
1970s, with the well-documented experiments at Stanford by 
Dr. John M.J. Madey (professor at that university until 
1988) . The EEL by itself has a background which dates back 
to at least 1947 and a patent by Elmer Corn, in which he 
describes amplification producing fast electromagnetic wave 
(periodic electron beam) interactions [2]. 
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In the 1950s, Motz worked on an undulator, and he 
passed a beam from a 3-MeV accelerator through it to produce 
millimeter wave radiation. He later used a 100-MeV beam to 
produce visible light. Another experiment of Motz was the 
Ubitron, a mildly relativistic microwave FEL. The power 
capability of this device was so much greater than that of 
conventional klystrons and magnetrons that records for peak 
power generation, which would stand for two decades, were 
established at both centimeter and millimeter wavelengths 
[ 2 ] . 

Nowadays the tendency of the institutions is to design 
EELS to operate at shorter wavelengths. The most powerful 
EEL in the United States is located in Jefferson's Lab at 
Newport News, Virginia; having a power of 14 kilowatts, this 
device works in an oscillator configuration. 

C. ADVANTAGES AND DISADVANTAGES OF THE FREE ELECTRON 

LASER 

The Free Electron laser has several advantages over 
other lasers. It is more flexible than the others; for 
example, it is tunable by design. That means that the output 
wavelength can be changed as desired by changing the energy 
of the electron beam or the properties of the undulator. A 
second advantage is that the EEL can operate continuously 
(24 hours a day) at a relatively low cost, whereas chemical 
lasers have to be replenished after a short time of 
operation, and they may generate toxic exhausts. The FEL 
will operate all the time if it has a continuous power 
supply and will not damage the environment with a toxic 
exhaust. Also, EEL can reach higher powers than the chemical 
lasers, because their lasing medium does not get hot. 
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II. BASIC COMPONENTS OF A FREE ELECTRON LASER 


There are two basic configurations of a free electron 
laser: The oscillator configuration and the amplifier 
configuration. Both configurations have many common 
components; the difference is in the way that the laser beam 
is generated. The basic components are: injector, linear 
accelerator, undulator, resonator (or, for the amplifier, a 
laser seed), beam dump and accessories. Figure 1 shows the 
basic components of a free electron laser. 


Electron Accelerator 



Figure 1. Basic components of a free electron laser 

(oscillator mode) from [3]. 


5 



















A. 


THE INJECTOR 


This is the first part of the device and is the source 
of the electron beam. Basically it contains a cathode, as 
does a TV tube. The electron beam is generated in short 
pulses (on the order of picoseconds) with energy of 5 MeV 
approximately, a value common to both configurations. There 
are several types of injectors; among these we can find the 
direct-current injector (DC), the radio-frequency injector 
(RF) and the semiconducting radio-frequency injector; and 
all of them can be used in an FEL system. Some injectors are 
driven by a laser in order to provide the required beam; the 
drive laser consists of a diode-pumped NdiGlass oscillator 
and a chirped pulse amplification system consisting of a 
granting stretcher, a flashlamp-pumped Nd:Glass regenerative 
amplifier and a grating compressor [8] . Figure 2 shows a 
basic scheme of a photocathode injector. 



Figure 2. 


Basic scheme of a photocathode injector from 

[4] . 
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B. LINEAR ACCELERATOR 

The linear accelerator (linac) is a critical component 
of the equipment, and increases the energy of the electron 
beam from 5 MeV at the output of the injector to about 100 
MeV. This increment of energy is done in the accelerator RF 
cavities. Electrons gain energy from the RF fields as they 
go from one cavity to the next, traveling through an 
electric field. The efficiency and gradient of the linac, 
and therefore the entire system, can be increased by use of 
a superconducting accelerator, which must operate at 
cryogenic temperatures. This requires the installation of a 
liquid helium refrigerator. Figure 3 shows the schematic of 
a linear accelerator. 





Figure 3. 


Schematic of a linear accelerator from [4]. 
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c. 


UNDULATOR 


The undulator is where the laser beam is created and 
enhanced, and where much of the interesting physics occur. 
An undulator is a series of alternating magnets, typically 
forming about 20 - 100 periods, separated by a few 

centimeters, and delivering a typical magnetic peak field of 
about 1 Tesla. The characteristics of the laser beam are 
determined by the undulator design, being a region in which 
the electrons are made to wiggle by a static alternating 
magnetic field. Since wiggling electrons radiate light, we 
can store this light in a resonator; and we can stimulate 
emission by radiating in the presence of this same light. 
While electrons are in the presence of the light there is an 
exchange of energy to and from the light. The electron 
distribution in phase space must be arranged in order to 
obtain the most efficient transfer of energy from the 
electrons to the light. Additional information will be 
written about this effect in a later chapter. In the 
amplifier configuration, the light passes only one time 
through the undulator, and then we require a longer 
undulator to extract the largest possible energy from the 
electron beam. Figure 4 shows the basic architecture of an 
undulator. 
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D. 


RESONATOR 


The resonator is the component where the light is 
stored in the oscillator configuration. It consists of an 
evacuated cavity bracketed by two mirrors typically 

separated by 20 - 30 meters, between which the light is 

stored. One of the mirrors is perfectly reflective (100%), 

the other is partially transmitting (it transmits some 
percentage of the light). It is through this mirror that the 
laser beam escapes to the outside. Resonators are often 
characterized by their quality factor Q, which is 
established by the inverse of their loss per pass (1 / Loss 
per pass). Figure 5 shows a basic arrangement of a 

resonator. 


electron beam pulses 
^ wiggler 




100 % 

reflecting 

mirror 


• •> 

• ^ 

partially 

transmitting 

mirror 


optical 

laser 

pulses 


Figure 5 . 


Schematic 


of a resonator from [3] . 


E. SEED LASER 

For an amplifier configuration it is necessary to have 
an external laser source, because unlike the oscillator, no 
light is stored in the optical cavity. This external laser 
source is the "seed laser"; it generates a laser beam of 
weak field (low power usually 10 - 100 W) that is introduced 
to the undulator for amplification to the desired output (1 
MW in our case) . Figure 6 shows a basic arrangement of a 
seed laser. 
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High-Power 
PEL Beam 



Figure 6. Arrangement of a seed laser from [5]. 

F. ACCESSORIES 

"Accessories" implies all the devices that control the 
behavior of the beams (electron and optical). They consist 
of bending and focusing magnets that center the electron 
beam and change its direction, and the mirrors that direct 
the optical beam. 

G. ELECTRON BEAM DUMP 

It is the last part of the electron path. In an FEL 
with recirculation of the electron beam, the electrons enter 
to the linac a second time 180° out of phase with the RF 
fields, giving most of their energy back to the electric 
field. At the output of the linac, the recirculated electron 
beam has energy of about 5 MeV, which is about the same 
energy that they have when they are generated at the 
injector. These electrons, with the residual energy, are 
dumped into what is called the "electron beam dump." At that 
relatively low energy, the electrons create less heat and 
radiation, so less shielding is required on this electron 
beam dump. 
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III. BASIC THEORY OF A FREE ELECTRON LASER 


A. THE PENDULUM EQUATION 

Since an accelerated charged particle emits radiation, 
we can create radiation by accelerating electrons inside the 
undulator. We do this by making the relativistic electrons 
'wiggle' through the undulator. When an electron enters to 
the undulator, it goes through electric and magnetic fields 
and encounters a force that acts upon it given by the 
Lorentz force equation. For relativistic electrons in the 


undulator, the 

complete 

Lorentz FEL force 

equations are 

given 

by [3] 

d[r~p) 

= -—(e+~pxb], 

me ' ' 

(III.l) 



dt 



dy 

dt 

= - — JXE, 
me 

(III .2) 




1 —2 

r 

(III.3) 

where 

e is the 

electron 

charge magnitude. 

m the electron 


mass, E the electric field, B the magnetic field, y0 = v/c. 


V is the electron velocity, c is the speed of light and / 
is the relativistic Lorentz factor. 

If we substitute the electric and magnetic fields of 
the FEL into equation (III.l), we get 


d(r/^±) 

dt 


— £’(l-y5j(cos^y,-sin^,0) + ;5^5(-sin(fc^z),cos(fc„z),0) (III. 4) 

yyj/^ L k / _ 


dt 


—(A cos y/ -smy/) + Bsm{k^z)-cos{k^z)) 

mr ^ ■ / \ / j 


(III.5) 
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where P^= is the magnitude of the longitudinal 

z component of the electron velocity, E is the magnitude of 
the electric field, B is the magnitude of the magnetic field 
and xj/ is the phase. We can solve for the transverse 

velocity components =-(^/ 7 )(cos(^^z),sin(^^z), 0 ^ where [5^ and 

Py are the components of the velocity and K = eB\llnmc^ . 

Substituting the transverse velocity components 
into (III.2) we get 

° e - - cKE cKE 

y = - P'E = -(cosA:^zcos^y-sinA:^zsin^) =-cos(4" + ^), (III. 6) 

me ymc ymc 

where P = {kQ + k^z-(ot is the electron phase, which represents 

the position of the electron within an optical wavelength. 
If we take the derivative of the electron phase with respect 
to the dimensionless time r = c?/L, where L = N\ is the length 

of the undulator, and N the number of periods of the 
undulator, we get the electron phase velocity 


v = ( = L[(k + k^)fi^-k]. 

For highly relativistic electrons (f»l), we can approximate 

P^^l-(K^+l)l2y\ and near resonance {k»kP^, so that the 
derivative of the electron phase becomes 


P = L [k + kP)(\-[K^+\)l2y^Y^ =Lk(\-[K^ + \)l2y^Y 


If we take the second derivative of the electron phase, we 
get 




yL^k 
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Substituting (III.6) into the last equation and defining 
\a\ = A7T:NeK\E\LlY^mc^ as the dimensionless field amplitude, we 
obtain the dimensionless pendulum equation. 


= K = |a|cos(i^ + ^z)) . (III.7) 

This equation describes the microscopic motion of electrons 
in phase space within the undulator. 

B. PHASE SPACE 

The phase space representation is a diagram that shows 
the evolution in position and velocity coordinates. It is 
useful to show, as in the last section, the electron energy 
exchange in an FEL. This diagram shows the electron phase 
velocity versus the phase. There are two types of orbits on 
a phase space diagram: closed and open orbits. 

Closed orbits for a classical pendulum represent a 
pendulum swinging back and forth; open orbits occur when the 
pendulum goes over the top due to a large angular velocity. 
The border between those two states is called the 
separatrix. Figure 7 shows a separatrix (black line) and the 
closed (green line) and open (red doted line) orbits on a 
phase space diagram; horizontal axis is the angular 
displacement and vertical axis the angular velocity. 


13 






Figure 7. Separatrix and orbits on a pendulum phase 

space diagram. 

Since we suggested the similarity between the FEL and a 
pendulum, and we obtained a pendulum equation describing the 
behavior of electrons in an FEL, we can state that \a\ and (j) 

are roughly constant for the case of low gain in the FEL; in 
the same way, the electron phase 4" can be viewed as the 
electron's position within an optical wavelength, and v is 
its phase velocity relative to resonance (v=0). 

For the FEL, the separatrix is given by [3]: 

= 2|a|[l + sin(4' + ^zJ)] (III.8) 

Figure 8 shows a phase space diagram of the behavior of 
approximately 20 electrons evenly spaced in ^ with vo=0 in a 
free electron laser. Yellow color indicates the electrons at 
the beginning of the undulator (t=0) and blue color 
indicates the electrons at the end of the undulator (t=1) . 
In this particular situation the gain of the system is zero, 
because half of the electrons gain energy and half of the 
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electrons lose it. In order to get gain, it is necessary to 
initiate the process with the electrons off resonance as 
Figure 9 shows. 


5 


V 


-5 


Figure 8. FEL phase space diagram from [3] 

C. ELECTRON BEAM AND OPTICAL BEAM INTERACTION 

One important issue in an FEL is the energy exchange 
between electrons and the optical field. As noted 
previously, values of v close to resonance (v=0) are 
required for optimal energy exchange. Phase space evolutions 
of the electrons showing a good bunching or grouping result 
in bigger gain. So it is necessary to create an optimal 
bunching to get the best values of gain. 

Figure 9 shows a phase space diagram of the evolution 
of the electrons in a free electron laser for a more 
realistic beam with 1000 particles and a small initial 
energy spread. 
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-12 

-31/2 J; 331/2 


Figure 9. FEL electron's bunching from [3]. 


D. FEL WAVE EQUATION 


From Maxwell's equations, we can obtain a full wave 
equation for the FEL optical field in terms of the 
transverse current [3]: 


V"A- 


1 An 


J, , 


(III.9) 


c" df c 

where A is the optical magnetic vector potential and J is 
the current density. 

For our helical undulator, a form for the optical 
vector potential can be written as [3]: 


A(x,?) = -!^£e'“ 


(III. 10) 


where E(z,t) is the complex optical amplitude, a = kz-cot is 
the carrier wave phase, i = (-/,!, 0) is the optical field's 
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polarization vector for circularly polarized light, and 
X = ri + yjf zk is the position vector. 

Substituting the solution in the wave equation we can 
get the results for the four dimensions. We can make some 
simplifications and substitutions like the slowly varying 
amplitude and phase approximation that will be discussed 
later in this chapter, in order to get the FEL wave equation 
as [ 3 ] : 

T 

dr 

where a(r) = |a|e'^ is the complex dimensionless optical field 
amplitude; the dimensionless current is j = %N{e7T:KLf'n^j y^mc^ , 
where is the beam particle density, measuring the 

interaction between the electron beam and the optical mode; 
and <e^‘^ > measures the amount of electron bunching. If j is 
large (J >> n) , the optical field a changes rapidly, so we 
have high gain. If j is small (J ^ n), we have low gain. 


)_ 


< -je"^ > 


(III.11) 



E. 


WEAK FIELD GAIN 


As mentioned previously, we can obtain gain through the 
interaction of the bunched electron beam with the optical 
field inside the undulator. An optical field is considered 
weak when |a|<;r it is smaller than n. The following wave 
equations are valid for weak and strong optical fields [3]: 

d\a\ 


dr 


= - j + (j))) , and 


(III.12) 


^ = y(sin(^ + ^)) . 


(Ill. 13) 
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They describe the change in optical amplitude \a\ and 

phase cp during FEL interaction with the electron beam. The 
pointing brackets denote the average of all the electrons 
sampled within an optical wavelength of the beam. 

From the coupling between the electron beam and optical 
field, we can calculate the dimensionless current density j. 


1. Low Gain Regime (j ^ n) 

In this regime, since j is small there is not an 
appreciable change in the optical field amplitude or phase 
as the electrons pass through the undulator. The evolution 
of the optical field amplitude and phase versus i 
(dimensionless time for the electrons to pass through the 
undulator from t=0 through t=1) can be obtained by expanding 
^ and V in powers of |a| = ao to get [7] 


(III.14) 


InlrW-n 

1 + / 

^2-2cos(v'or)-KQrsin(v'or)^ 


J 

1 




2sin (kot) - VqT (l + cos (vor)) 


(Ill.15) 


where is the initial optical field amplitude, and is 
the initial phase velocity of the electron beam. 

In an FEL system, the gain is sensitive to the initial 
phase velocity Vq . If the electron beam starts with ^0=0 
(at resonance), the resultant gain is negligible. The 
resulting gain inside the undulator assuming weak fields is 
given by [7] 

2 - 2cos (Kgr) - sin (i/gr) ^ 


G(r) = j 


Vn 


(III.16) 
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2. High Gain Regime (j»n) 

In this regime, the dimensionless current density is 
large, and it strongly affects the optical field amplitude 
and phase as the electrons go through the undulator. In this 
case, the optical field amplitude and phase during the 
electron beam interaction with the optical field are given 
by [7] 

f 

|a(r)| = —^ r (III.17) 

I V ;i 3 



(III. 18) 


From these formulas, we can see large changes in the 
amplitude and phase due to the current density; then for 
large changes in the phase, the separatrix is shifted and 
the optimum bunching location changes on the phase space 
diagram. In contrast to the low gain regime, the gain is not 
as sensitive to the initial phase velocity, but the gain 
spectrum changes as well. 

The expression for the gain in the undulator for this 
regime is given by [7] 



(III. 19) 


F. STRONG FIELD GAIN 

For this case, in both regimes (low and high gain) , 
strong optical fields make the electron beam overbunch in 
the undulator as the electrons travel in it. This process is 
called saturation. Overbunching creates absorption of energy 
by the electrons from the optical field, reducing the gain 
of the system. In the low gain regime, saturation is 
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achieved at much lower initial field amplitudes than 

the high gain regime [V]. 

When there is not saturation, the phase space evolution 
of the electrons is slower and the bunch of electrons travel 
downward to a phase of ^ ~ n. However, the stronger optical 
field causes a greater force on the electrons and more 
evolution in phase space. The bunch forms and continues 
along closed orbits until it begins an upward trend, taking 
the energy away from the optical field. 

G. TAPERED WIGGLER 

In this work, an analysis of an FEL in the amplifier 
configuration is made; this type of system requires a much 
longer undulator than the oscillator. If we want to increase 
the gain and extraction of the system, we can extend the 
saturation of electrons by tapering the undulator. Taper is 

a change in the resonance condition X = + as the 

electrons travel along the undulator; this alteration 
extends the energy exchange between electrons and the 
optical field; this extended exchange allows a greater 
amount of energy taken from the electrons and an increased 
extraction of energy from the device. As the electrons 
evolve, they loose energy, and y decreases, then, in order 
to keep the same wavelength, or K must be decreased. 

There are several ways to taper the undulator; the most 
common is to vary the undulator parameter K by changing the 
magnetic fields along it. The easiest way to do this taper 
is to change the size of the gap between magnets by a small 
amount (Tenths of millimeters) ; this will change the 
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strength of the field at certain points. In this way, we can 
generate different types of tapering. 

1. Linear Tapering 

This tapering method consists of a continuous increment 
or decrement of the magnetic field strength along the 
undulator. This variation of the field can be negative 
(increasing gap) or positive (decreasing gap). The taper can 
start at any location along the undulator and continues 
linearly through the rest of the undulator. Typically, the 
taper stars near the location where the optical field is 
expected to reach its saturation value . 

As we can see in Figures 10 and 11 (positive and 
negative linear taper), this technique exposes the electron 
beam to a linearly changing magnetic field. 



Figure 10. Positive linear tapered undulator from [3]. 
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Figure 11. Negative linear tapered undulator from [3]. 

2. Step-Tapering 

Other useful taper method is the step taper. This 
method consists of changing the gap at a certain point 
within the undulator and maintaining this gap through the 
rest of it. Figures 12 and 13 show a positive and a negative 
step-tapered undulator. 



Figure 13. 


Negative step-tapered undulator from [3]. 
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As the electron beam's energy changes along the 
undulator, the resonance condition is altered. By tapering 
the undulator, we introduce the following phase acceleration 
in the pendulum equation [7]: 

, , AK 

S = -47tN -^ (III.20) 

l + K^ K 

where AK/K is the fractional change on the undulator 
parameter, and where K is proportional to the magnetic 
field. 

The pendulum equation for the tapered undulator becomes 

[7] : 


= v = S+a cos[<^ + (p) . 


(Ill .21) 


The optimal extraction from a normal undulator is [7]: 




1 

2N 


but in a tapered undulator, the phase acceleration 
an important role in the extraction. The extraction 
case is [7] : 


(III.22) 

5 plays 
for this 


1^1 

n. « -L-L . (111. 2 3 ) 

^ SttN 

As the phase acceleration (positive or negative) 
approaches the value of the optical field amplitude, the 
separatrix becomes smaller. When this acceleration is 
greater, the separatrix vanishes, and bunching does not 
occur because the electrons are not trapped. 

It is obvious that at bigger phase acceleration, the 
extraction will be bigger; but there is an upper limit. This 
limit is set by the fact that the pendulum equation has no 
solution when the taper phase acceleration is bigger than 
the optical field amplitude; however, a lower limit is 
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determined by the maximum deceleration a tapered electron 
can undergo in the no tapered case. Then, we can state the 
taper range as: 

4Va < <5’ < |a| 

Figures 14 and 15 show the phase space evolution of a 
positive and negative linearly tapered undulator. In the 
positive taper we can see that some of the electrons are 
trapped inside the separatrix and they bunch, improving the 
energy exchange. In the negative taper case, we can see that 
electrons go around the separatrix, but we still have 
bunching and good extraction. For this last case, the taper 
will work well only when electron beam starts above 
resonance. Optimal taper occurs when approximately half of 
the electrons are trapped. 



Figure 14. 


Phase space diagram of a positive linearly 
tapered undulator from [3]. 
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Figure 15. Phase space diagram of a negative linearly 

tapered undulator from [3]. 

Figures 16 and 17 show the phase space evolution of a 
positive and negative step-tapered undulator. There we can 
see that positive taper traps electrons like the positive 
linearly tapered case; and negative tapered electrons are 
similar to the negative linearly tapered case. 



Figure 16. 


Phase space diagram of a positive step- 
tapered undulator from [3]. 
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Figure 17. 


Phase space diagram of a negative step- 
tapered undulator from [3]. 


H. DIFFRACTION 


Propagation of a laser beam can be described in free 
space by Maxwell's wave equation [7] 


C 


2 a^2 

C Ot 


V 


a(x,t) = 0 , 


(III .24) 


where a{x,t)e'“ = 


a(x,t) 


i^a{x,t) ia 


e is the optical electric field. 


For laser light, the beam is coherent and the electric 
field is slow-varying in the direction of propagation over 
an optical wavelength and slow-varying in time over the 
optical frequency. In the last equation, a= kz-cot represents 
the phase of the carrier wave, and the optical field is 
represented by a, as defined earlier. 

If the optical field is constant, the electrical field 
describes a simple plane wave propagating in the z 
direction, containing only one frequency co. If the wave's 
amplitude or phase contain spatial and time dependence the 
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laser beam will have a finite-with spectrum; if both 
dependences are slow then the laser will have a narrow 
spectrum. In an FEL, the amplitude and phase are slowly- 
varying, so we can make the following assumptions: 

d\a\ , , d(A , d\a\ , , dd 

, —<^o)(b , —^«,k\a\ , —<^k(h (III. 2 5) 

dt 5? dz 

After substitutions we get 


f d^a 5a 2 1 -2 

Via+ — :r + 2ik - ka -c 

dz^ dz 


d^a da 2 ^ ia r. 

— --2ico - 0 ) a e =0 

dt^ dt 


(III.26) 


Using assumptions (III.25), we can neglect the terms 

with two derivatives. Multiplying by and using co=kc we 

arrive at 

f d 15^ ^ 

V ±+2ik 1- a{x,t) = 0 . (III. 27) 

5z c dt) 

This equation can be simplified further by introducing 
a coordinate transformation u=z-ct and a dimensionless time 
T=ct/L where L is the length of the undulator. Then, the 
wave equation can be written with one less partial 
derivative, 

a(x,t) = 0 , or (III.28) 


2ik d 

Vi +- 

L dr 


-- a(x,t) = 0. 

2k dr 


(III .29) 


This equation is known as the parabolic wave equation. 

The first term in this equation is the diffraction term, and 

its magnitude depends on the transverse area of the laser 

beam, the wavelength (X=2n/k), and the undulator length L. 

The Rayleigh length (distance from the mode waist to 

where the area of the optical mode doubles due to 

diffraction) relates the optical mode waist radius to 
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the wavelength of the light by Z^-TiWllX. In dimensionless 
units Zo=Z^/L. Based on diffraction along the undulator, so 
the characteristic mode waist radius is 
■\/ LX 17t = -sIlLI k . If the optical mode radius is much larger 

than \IlZITT , then the diffraction has a small effect over 
the undulator. If it is smaller, diffraction has a big 
effect on propagation, changing both the amplitude and phase 
of the optical field. 

It is important to mention that this description of 
diffraction is calculated for laser light in free space. For 
an FEL in the amplifier configuration, diffraction is 
altered by the electron interaction. Inside the undulator, 
the wave equation is inhomogeneous, and the Gaussian 
distribution of the optical field changes; the laser beam is 
affected by the electrons that try to focus the laser beam. 
Then we have to consider the last term of the FEL wave 
equation <-je'^>. This new term affects the diffraction of 
the optical field. 

I. SHORT PULSES AND SLIPPAGE 

Eor energy exchange between the electrons and optical 
field, they must overlap inside the undulator. During their 
pass through the wiggler, electrons and photons have 
different velocities. Photons travel at speed of light, and 
relativistic electrons move close to it but slower. The term 
used to describe the overlap of the optical and electron 
pulses is known as "slippage"; according to the resonance 
condition, the electrons drop behind the light by a slippage 
distance NX over the length of the undulator. Eigure 18 
shows the slippage concept. Red pulse is the electron pulse 
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and blue pulse is the optical pulse; the left overlap 
represents the position at the beginning of the undulator, 
and the right overlap is at the end of the undulator. 



Figure 18. Electron and light slippage from [3] . 


We have short pulses when the micropulse length is 

on the order of or less than the slippage distance NX. In 
order to get the required energy transfer, the short 
electron pulses and optical pulses must be synchronized. 
That process develops in stages, and electron bunching must 
take place in the right moment. At the beginning of the 
undulator the electrons and light overlap, but there is no 
bunching and therefore no gain. At the middle, electrons 
during their bunching process slightly amplify the light 
pulse. At the end, bunched electrons amplify the trailing 
part of the optical pulse. As a result, the optical pulse 
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centroid travels slower than speed of light; this process is 
called "optical lethargy". Figure 19 shows this process. 



Figure 19. Optical lethargy from [3]. 
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IV. SIMULATIONS OF FREE ELECTRON LASERS IN 
AMPLIFIER CONFIGURATION 


In this chapter, several free electron laser amplifiers 
were explored with computer simulations, using the NFS Apple 
Xserve cluster with 64 nodes/128 processors. On this device, 
a simulation program was run that solved the coupled 
pendulum and FEL wave equations in order to understand the 
four dimensional behavior of the high power FEL amplifier. 
The simulation program required dimensionless input 
parameters, which make the results applicable to a broad 
range of FELs. 

A. DESCRIPTION OF THE INITIAL PARAMETERS 

The first step in this process is to establish the 
initial parameters that describe the system to be evaluated. 
Using these parameters, the dimensionless input variables 
are calculated in order to run the simulation. 

We can classify these parameters as primary and 
secondary. Secondary parameters are coupled to the primary 
parameters and they have to be calculated using several 
formulas that will be described in this work. 

The desired characteristics of the device were 
introduced into an Excel spread sheet, which calculated the 
secondary parameters using formulas, and the dimensionless 
input parameters for the simulation program. Tables 1, 2 and 
3 show the description of the dimensional parameters that 
were used and the formulae used to obtain the secondary 
dimensional parameters. 
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■undulator parameters 


\) 

Undulator period 


N 

Number of periods 


g 

Undulator gap 


B 

Magnetic field (rms) 


K 

Undulator parameter 

iv — /_ 2 

/ 2nmc 

L 

Undulator length 

< 

II 


Table 1. Undulator parameters used to describe an FEL in a 

simulation. 


ELECTRON BEAM PARAMETERS 


<lb 


Q 


r 

^ Peak 
^Avg 

A/ 

r 

AO 


Beam kinetic energy 
Bunch charge 

Beam radius (rms) 

Pulse duration, FWHM 

Pulse repetition 
frequency 

Pulse length, FWHM 
Lorentz factor 

Peak current 

Average current 
Normalized rms emittance 
Longitudinal emittance 

Beam energy spread (rms) 

Beam angular spread 
(rms) 

Beam average power 




'Ay 


4 4^ 


_(£,+0.511Mey) 


= % 
Peak / f 



^Avg ~ 


^ ^ 235£i/ 
y /(^b)(h) 

A6» = ^«/ 

/rh 

P = F I 

^b ^ Avg 



Table 2. Electron beam parameters used to describe an FEL 

in a simulation. 
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OPTICAL PARAMETERS 


Seed laser power 


T 

in 

S 




Seed pulse duration 


Distance to First Optic 


Rayleigh length 


X Optical wavelength 


Z 

^0 “ /Z 


ly^ 



Mode waist radius, 1/e 





Table 3. Optical parameters used to describe an FEL in a 

simulation. 


B. SIMULATIONS OF THE FIRST FEL 

This first laser was intended to have a power of 100 kW 
with a wavelength of 1.6 pm in an amplifier configuration. 
The undulator has 120 periods, and its length is 324 cm. 
Tables 4 through 7 show these parameters. 


I ELECTRON BEAM PARAMETERS 

E, 

Beam energy (MeV) 

80 


Bunch charge (nC) 

0.2 

h 

Beam radius, (mm) 

0.13 

h 

Pulse duration, FWHM (ps) 

1.0 

Q 

Pulse rep frequency (MHz) 

703 

4 

Pulse length, FWHM (cm) 

0.03 

Y 

Lorentz factor 

158 

^ Peak 

Peak current (A) 

200 

^Avg 

Average current (mA) 

141 


Normalized rms emittance (mm mrad) 

3.5 


Longitudinal emittance (keV ps) 

36 

<] 

Beam energy spread (%) 

0.11 

A0 

Beam angular spread (mrad) 

0.17 

n 

Beam average power (MW) 

11 


Initial electron beam parameters of the first 
simulated FEL device. 
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Table 4. 

























UNDULATOR PARAMETERS 


Undulator period (cm) 2.7 
N Number of periods 120 
8 Undulator gap (cm) 1.00 
K Undulator parameter, rms 1.39 
L Undulator length (cm) 324 


Table 5. Initial undulator parameters of the first 

simulated FEL device. 



Table 6. Initial optical parameters of the first simulated 

FEL device. 



Table 7. Dimensionless parameters of the first simulated 

FEL device. 
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The first parameter to analyze was the taper start time 
. This time is established as a dimensionless parameter 
using a scale from 0 to 1, with ^"^=0 corresponding to the 
beginning of the undulator and ^^=1 the end of he 
undulator. The taper start time was simulated for several 
values from 0.4 to 0.8 with intervals of 0.1; using a taper 
rate of 5=40n. The start time with the greatest extraction 
was at r, = 0.7 with an extraction g = 0.92%. Figure 20 shows 
a graph of the obtained extractions versus the simulated 
taper start times for 5=40n. 



Figure 20. Graph of extraction vs. taper start times. 

The full simulation output for the most efficient taper 
start time (= 0.7 ) at 5=40n is shown on Figure 21. There we 
can see the power P(t) growing to saturation at the end of 
the undulator (t=1), with about half the electrons trapped 
in phase space (^,v) as shown on phase space diagram. Also 
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we can see a typical diffraction pattern of the beam in x 
and y. Slippage produces a narrow optical pulse in the z 
direction. 


Simulation run time: 39.5 sec Fri May 4 11:31:22 2007 

nx=200, nz=33, nt=2 5G, n|i=100000, Wp=12, seed=7, wbins=33, ebins=16 



2 -2 


2 -4 


4-13 


13 



0.0 X 1.0 0m7-2x 2-5x 5 

F=0.44, G=1.14e+03, ti =0.00921, ^7/7=0.0408 
THEORY: G(jF»l)=5. 14e+04 
THEOm: r|(jF»l)=0. 0125 

a =141 at X =0.70, THE»RY: a =164 at x =0.57 
s s s s 

M^=0.047, c^(0,0)=0.25, (0,2)=0.11, c^(2,0)=0.11 


Figure 21. 


Full simulation output for taper start time 

0.7. 
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Once the most efficient value of the taper start time 
(r^=0.7) was found, it was used to run a set of simulations 
to analyze the dimensionless linear taper rate (5) . 
Simulations were run for 5 values from 0 to 80n at intervals 
of lOn. The most efficient rate was at 5 = 40n with an 

extraction r] = 0.92%. Figure 22 shows a graph of the 

obtained extractions versus simulated taper rates. 



Figure 22. Graph of extraction vs. linear taper rates. 


The full simulation output for the most efficient taper 
start time (rj=0.7) is shown in Figure 21 (coincidently the 
taper rate used in that simulation was the most efficient). 
The results are the same as the previous simulation. Since 
the extraction is not enough to achieve 100 kW, it is 
necessary to change other parameters in order to reach the 
goal. 
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In order to look for a better extraction, a longer 
undulator was tried, increasing the number of periods to 
140. Since the length of the undulator changed to 378 cm, it 
was necessary to calculate the dimensionless input 
parameters again. Table 8 shows only the parameters that 
changed for this case. 


I ELECTRON BEAM PARAMETERS 

A0 

Beam 

Beam 

radius, 
angular 

(mm) 

spread (mrad) 

0.14 

0.15 







Table 8. New parameters for N=140. 


Using a taper start time r^=0.5, a set of simulations 
explored the effect of the dimensionless taper rate (5). 
Simulations were run for 5 values from 0 to 80n at intervals 
of lOn. The most efficient rate was at 5 = 30n with an 
extraction g = 1.13%. Figure 23 shows a graph of obtained 
extractions versus the simulated taper rates. 
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Dimensionless Linear Taper Rate in Multiples of pi 

Figure 23. Graph of extraction vs. linear taper rates 

for N=140. 


With this longer undulator we can see that the 
extraction is bigger even with a smaller taper rate. 
Diffraction patterns are very similar to the previous case 
and a higher power reaches saturation at the end of the 
undulator. On Figure 24 we can see that electron bunching is 
quite similar to the previous case, but the optical field 
evolved into a combination of higher-order modes. 



-11/2 ^ 3ii/2 -2x2 


Figure 24. Electron bunching and optical field for 

linear taper rate 5=30n. 
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Continuing with this exploration, and using the optimal 
taper rate found for this FEL (5=30n) and a taper start time 
t^=0.6, this device was analyzed by changing the electron 
beam pulse duration () in a range from 0.4 to 1.6 ps in 
0.2 ps steps. It is important to mention that any change in 
this parameter changes the input parameters of the 
simulation, so a recalculation of dimensionless parameters 
was made for every value of . The largest extraction (g = 
1.18%) was found at t^=0.Sps. Figure 25 shows a graph of the 
obtained extractions versus the simulated electron beam 
pulse duration. 



Figure 25. Graph of extraction vs. electron beam pulse 

duration. 
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The results at the best pulse duration (?^=0.8ps) are 
similar to the last case, but this time it did not have many 
higher-order optical modes. 

The next step is to change other parameters to obtain a 
better extraction. This time the bunch charge was changed 
from 0.2 to 0.4 nC, with the pulse repetition frequency cut 
in half (Q = 350 MHz) to keep the average current the same 
(140 mA) . This change affects some of the input parameters 
to the simulation program, so, it is necessary to 
recalculate the dimensionless parameters to create another 
input file. Table 9 shows the initial parameters that 
changed for this new case with respect to the last case upon 
changing the bunch charge value. 




Table 9. Parameters that changed for micropulse charge of 

0.4 
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Using these new parameters the input file was updated 
for the simulation, and several runs tried to obtain optimal 
outputs for the same ranges of taper values. 

Again the first parameter to analyze was the taper 
start time. The taper start time was simulated for several 
values from 0.4 to 0.8 at intervals of 0.1; using a linear 
taper rate of 40n. The optimum start time was at 0.5 with an 
extraction g = 2.06%. Figure 27 shows a graph of the 
obtained extractions versus the simulated taper times. 



Figure 26. Extraction vs. dimensionless taper start time 

for micropulse charge of 0.4 nC. 

The results for the optimum taper time (0.5) are very 
similar to the last case, but this there was slightly better 
electron bunching, resulting in a better extraction. As a 
result and a higher power was delivered when saturation was 
reached at the end of the undulator. 
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Again, once the optimum value of the taper start time 
(0.5) was found, it was used to run a set of simulations to 
analyze the dimensionless taper rate (5) . Simulations were 
run for 5 values from 0 to 80n at intervals of lOn. The 
optimum rate was at 5 = 60n and 5 = 70n with an extraction p 
= 2.32%. Figure 29 shows a graph of the obtained extractions 
versus the simulated taper rates. 



Figure 27. Extraction vs. dimensionless taper rates for 

micropulse charge of 0.4 nC. 

The results for the taper rate 5 = 60n are shown in 
Figure 28. Extraction is much better in this case, but the 
taper rate must be large in order to reach this level of 
extraction. As in the last case, a higher power was 
delivered when saturation was reached at the end of the 
undulator. It may be possible to reach these levels changing 
other parameters. 
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simulation run time: 40.7 sec Tue May 8 13:41:02 2007 

rDC=200, nz=33, nt=256, rp=100000, Wp=12, seed=7, vftiins=33, ebins=16 


***** pel amplifier 

4D SIMULATION ***** L 

|a| 

max 

j=9137 0^=0.1 

ay=0.1 

0^=0.42 

a^=0.42 0^=1. 3 

v^=6.19 

°V7~^ 

tOp=7.8 Xp=0.5 

Ay^=0 

II 

o 

5C=0 N=140 

D=0 

5=G0x 

5 =0ii X =0 . 6 

a =9.7 

Q 

II 

t =0 z =0.09 

X =3 

o 

II 

o 

s s 

o 

a 

w o 

m 

m 


1 

X 

0 


-2 X 2-2 Y 2-2 z 2-25 v 25 




0.0 X 1.0 Om7-2x 2-5x 5 

F=0.44, G=:2.7e+03, ti=0.0232, A7/7=0.078 5 
THEORY: G(jr»l)=l. 23e+04 
THEORY: ti( jF»l)=0. 0214 

a =329 at x =0.60, THEORY: a =321 at x =0.42 
s s s s 

M^=0.025, c^(0,0)=0.24, (0,2)=0.11, c^(2,0)=0.11 


Figure 28. Results for micropulse charge of 0.4 nC and 

taper rate 5=60n. 
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The last explored parameter was the emittance; this 
parameter is mainly determined by the injector. Emittance 
changes have a direct effect on the initial parameters; and 
any change of this parameter needs to be applied to the 
spread sheet in order to calculate the new input parameters. 

For this simulation, emittance was adjusted to 5.1 mm 
mrad, and Table 10 shows the initial parameters that changed 
for this new case with respect to the last case upon 
changing emittance. 



Normalized rms emittance (mm mrad) 5.1 

lS.6 Beam angular spread (mrad) 0.23 


1 OPTICAL PARAMETERS 


VKg Mode waist radius, 1/e (mm) 

0.40 



Table 10. Parameters that changed after changing the 

emittance. 


With these parameters and using a taper start time 
r^=0.5, it was run a set of simulations to explore the 
dimensionless taper rate (5) . Simulations were run for 5 
values from 30n to 80n at intervals of lOn. The optimum 
taper rate was at 5 = 60n with an extraction g = 2.15%. 
Figure 29 shows a graph of obtained extractions versus the 
simulated taper rates. 
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2.5 



Dimensionless Linear Taper Rate in Multiples of pi 


Figure 29. Extraction vs. dimensionless taper rates for 

emittance=5.1mm mrad. 

This graph shows that with this value of emittance, 
taper rate variation over this range does not have a big 

effect on the extraction. For 5 = 60n, the results of the 
simulation are very similar to the last case. As in the 

last case, a higher power was delivered when saturation was 
reached at the end of the undulator. 

C. SIMULATIONS OF THE SECOND FEL 

This second set of simulations explores an FEL 
amplifier with a power of 1 megawatt and an optical 

wavelength of 1.6 microns. The same procedure is used for 

this laser, but this time the number of periods of the 
undulator was analyzed for several values in order to find 
the optimum. The general input parameters of this device are 
shown on Tables 11, 12, 13 and 14. 
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ELECTRON BEAM PARAMETERS 


E, 

Beam energy (MeV) 

75.5 

‘Ih 

Bunch charge (nC) 

1.5 

h. 

Beam radius, x (mm) 

0.31 

% 

Beam radius, y (mm) 

0.18 

h 

Pulse duration, FWHM (ps) 

1.1 

n 

Pulse rep frequency (MHz) 

703 

h 

Pulse length, FWHM (cm) 

0.034 

r 

Lorentz factor 

149 

^ Peak 

Peak current (A) 

1333 

^ Avg 

Average current (mA) 

1055 


Normalized rms emittance, x (mm 

mrad) 15.0 

^ny 

Normalized rms emittance, y (mm 

mrad) 5.0 


Longitudinal emittance (keV ps) 

70 

Af 

Y 

Beam energy spread (%) 

0.19 


Beam angular spread, x (mrad) 

0.32 

A^v 

Beam angular spread, y (mrad) 

0.18 


Beam average power (MW) 

80 


Table 11. Initial electron beam parameters of the second 

simulated FEL. 


UNDULATOR PARAMETERS 


/l„ Undulator period (cm) 

N Number of periods 

g Undulator gap (cm) 

K Undulator parameter, rms 

L Undulator length (cm) 


2.62 
150 
1.00 
1.30 
393 


Table 12. Initial undulator parameters of the second 

simulated FEL device. 
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OPTICAL PARAMETERS 


Z’,, Seed laser power (W) 100 

T.^ Seed pulse duration, FWHM (ps) 2.3 

S Distance to First Optic (cm) 1000 

Zg Rayleigh length (cm) 105 

A Optical wavelength (microns) 1.59 

Wg Mode waist radius, 1/e (mm) 0.73 

Table 13. Initial optical parameters of the second simulated 

FEL device. 


DIMENSIONLESS PARAMETERS 


j Normalized current density, linear undulator 20288 
cr^ Normalized beam radius, x 0.22 

(Jy Normalized beam radius, y 0.13 

Normalized beam angular spread, x 0.8 9 

Normalized beam angular spread, y 0.52 

Phase velocity spread due to energy spread 3.7 
cr^ Normalized pulse length 1.4 

ag Normalized initial optical field amplitude 4.3 

cr^ Normalized initial optical pulse length 2.8 

Zg Normalized Rayleigh length 0.27 

Wg Normalized mode waist radius 0.52 

co^ Betatron oscillation frequency 8.2 


Table 14. Dimensionless parameters of the second simulated 

FEL device. 
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For the first simulation, it was necessary to update 
the dimensionless input file for each value of N. The number 
of periods was analyzed from N=100 through N=200 at 
intervals of 10, obtaining the highest extraction g = 2.21% 
at N=140. Figure 30 shows a graph of the extraction versus 
the number of periods of the undulator. 



Figure 30. Extraction vs. number of periods of the 

second FEL. 

Figure 31 shows the results of the simulation for this 
set of parameters with N=140, which is the optimum value.we 
can see the power P(t) growing to saturation at the end of 
the undulator (t=1) with an intermediate bump corresponding 
to a synchrotron oscillation, and about half the electrons 
trapped in phase space (^,v) as shown on phase space 
diagram. 
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Simulation run time: 39.9 sec Tue Aug 7 13:05:59 2007 

nx=200, nz=33, nt=2 56, r5D=100000, Wp=12, seed=7, ’wbins=33, ebins=20 



-3 X 


3-3 Y 


3 -4 


4-13 


13 



F=0.42, G=:2.03e+04, ii=0.0221, A7/ic0.09 54 
THE0R5f: G(jF»l)=2. 65e+04 
THEOrar: Ti(jF»l)=0. 0179 

a =535 at x =0.50, THEJORY: a =371 at x =0.48 
s s s s 

M^=0.058, c^(0,0)=0.33, (2,0)=0.13, c^(0,2)=0.12 


Figure 31. Results for N=140 of the second FEL. 
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Using the undulator with N=140, the taper start time 
was simulated for several values from 0.0 to 0.9 at 
intervals of 0.1 with a taper rate of 40n. The optimum start 
time (greatest extraction) was at 0.5 with an extraction g = 
2.21%. Figure 32 shows a graph of the obtained extractions 
versus the simulated taper times. 



Figure 32. Extraction vs. dimensionless taper start time 

of the second FEL. 

The results for the best taper time (0.5) are shown in 
Figure 31. 

The taper start time of 0.5 was used to run a set of 
simulations to analyze the dimensionless taper rate (5) . 
Simulations were run for 5 values from 0 to 150n at 
intervals of lOn. The optimum rate was at 5 = lOOn with an 
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Figure 33 shows a graph of the 


extraction g = 3.55%. 

obtained extractions versus the taper rates. 



Figure 33. Extraction vs. dimensionless linear taper 

rates of the second FEL. 

The results for the optimum taper rate 5 = lOOn are 
shown in Figure 34. Extraction is much better in this case. 
As in the last case, a higher power was delivered when 
saturation was reached at the end of the undulator. Less 
than half of the electrons bunched inside of the separatrix 
as shown on the phase space diagram. 
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simulation run time: 39.3 sec Ned ftug 8 12:23:38 2007 

nx=200, nz=33, nt=256, i^=100000, Wp=12, seed=7, vA>ins=33, ebins=20 


***** fel amplifier 4d simulation ***** 


|a I max 


j=11905 

Ox=0.22 

Oy=0.l3 

a. =0.89 

cy^=0.52 0^=1. 5 
5^=0 N=140 

V =4.54 

a =3.4 

V7 

5=100ji 

8 

II 

•m] 

Xp=0.5 

AY^=0 

II 

o 

D=0 

5 =0k 
s 

X =0.5 
s 

a =3.9 
o 

a =3 
a 

t =0 z =0.27 

w o 

X =2.7 
m 

0 =0 
m 


a(x, 0,0)1.08e+03 

a(0,Y,0X.08e+03 

a (0,0, zj. . 08e+03 

A 

P(0,0,^^^^^ 




1 ; 

1 : 

I i 

II 




LI 


i 



a(x,0,0) 

3.89 

a(0,Y,0) 3.89 

a(0,0,z) 3.89 

P(Q,0,V). ^ . 


-3 


3 -3 


3 -4 


4-13 


13 


2 58 


-258 



P(x) 

4.52x10^ 

G(x) 



(;,v) 




. 


a(x,Y,0) 
x=l 


-x/2 t, 37t/2 -3 




0.0 


1.0 0 m 7 -3 


F=0.42, G=3.2e+04, ti=0.0355, A7/7=0.139 
THEORY: G(jF»l)=6. Ole+04 
THEORY: r| (jF»l)=0. 0446 

a =489 at x =0.50, THEORY: a =371 at x =0.48 
s s s s 

M^=0.094, c^(0,0)=0.35, c^(2,0)=0.14, c^(0,2)=0.12 


Figure 34. Results for 5 = lOOn of the second FEL. 
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D. 


SIMULATIONS OF THE THIRD FEL 


This third set of simulations explores an FEL amplifier 
with a power of 3 megawatts and an optical wavelength of 1.6 
microns. In this case, more realistic electron beam 
parameters were generated by an external computer program 
(Parmela). The number of periods of the undulator was 
adjusted in order to obtain the most extraction. The general 
input parameters of this device are shown on Tables 15 
through 18. 


ELECTRON BEAM PARAMETERS 


Beam energy (MeV) 75.5 

qf, Bunch charge (nC) 1.5 

Beam radius, x (mm) 0.2 8 

r^y Beam radius, y (mm) 0.18 

Pulse duration, FWHM (ps) 0.88 

Q Pulse rep frequency (MHz) 703 

4 Pulse length, FWHM (cm) 0.026 

7 Lorentz factor 149 

^Peak Peak current (A) 17 05 

Average current (mA) 1055 

Normalized rms emittance, x (mm mrad) 14.0 
£„y Normalized rms emittance, y (mm mrad) 5.6 

s^ Longitudinal emittance (keV ps) 60 

Af 

- Beam energy spread (%) 0.21 

r 

Beam angular spread, x (mrad) 0.33 

^Oy Beam angular spread, y (mrad) 0.21 

Pfj Beam average power (MW) 80 


Table 15. Initial electron beam parameters of the third 

simulated FEL. 
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UNDULATOR PARAMETERS 


Aq 

Undulator 

period (cm) 

2.62 

N 

Number of 

periods 

130 

g 

Undulator 

gap (cm) 

1.00 

K 

Undulator 

parameter, rms 

1.30 

L 

Undulator 

length (cm) 

341 


Table 16. Initial undulator parameters of the third 

simulated FEL device. 



Table 17. Initial optical parameters of the third simulated 

FEL device. 



Table 18. Dimensionless parameters of the third simulated 

FEL device. 
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For the first simulation, it was necessary to update 
the dimensionless input file each time the number of periods 
was changed. The number of periods was varied from N=100 
through N=200 in intervals of 10, obtaining the highest 
extraction of ri=3.15% at N=130. Figure 35 shows a graph of 
the extraction versus the number of periods of the 
undulator. We also considered rotating the undulator 90° by 
switching the x and y components of the normalized beam 
radius and angular spread. This change did not make any 
significant difference in the results. 



Figure 35. Extraction vs. number of periods of the third 

FEL. 
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Figure 38 shows the simulation results for this set of 
parameters with N=130, which is the optimal value. We can 
see the power P(t) growing to saturation at the end of the 
undulator (t=1) with an intermediate synchrotron 
oscillation, and about half of the electrons bunched in 
phase space (^,v) as shown on phase space diagram. 

Using the optimum undulator (N=130), the taper start 
time was simulated for several values from 0.2 to 0.8 at 
intervals of 0.1 using a taper rate of 80n. The optimal 
start time (greatest extraction) was at 0.5 with an 
extraction g = 3.14%. Figure 36 shows a graph of the 
extractions versus the taper start times. 



Figure 36. Extraction vs. dimensionless taper start time 

of the third FEL. 
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The results for the optimal taper start time (0.5) are 
shown in Figure 38. 

Again, once the optimal value of the taper start time 
(0.5) was found, it was used to run a set of simulations 
varying the dimensionless linear taper rate (5). Simulations 
were run for 5 values from 20n to llOn at intervals of lOn. 
The most efficient rates were at 5 = 80n and 5 = 90n with an 
extraction g = 3.14%. Figure 37 shows a graph of the 
obtained extractions versus the taper rates. Beyond the peak 
it is roughly flat. 



Figure 37. Extraction vs. dimensionless taper rates of 

the third FEL. 

The results for the taper rate 5 = 80n are shown in 
Eigure 38. 
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simulation run time: 31.9 sec Tue Aug 21 10:09:38 2007 

nx=300, nz=33, nt=2 5G, i5)=30000, ^3=16, seed=7, vAiins=38, eblns=15 


* * * * * 

FEL AMPLIFIER 

4D SIMULATION ***** L 

Ia| 

max 

j=1283G 

0^=0.14 

o =0.22 

Ofl^,=0.8G a,=1.3 

Vq=0.75G 
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V7 

c«p=7.1 

Xp=0.5 

o 

II 

Ae =0 5C=0 M=130 

yo ^ 

D=0 

5=80x 

5 =0ii 

o 

II 

a =4 

a =2.6 t =0 z =0.26 

X =2.5 

II 

o 

s 

s 

o 

a w o 

m 

m 


1 

X 

0 


-4 X 4-4 Y 4-3 z 3-17 v 17 


a(x, o,o; 

1.07e+03 

a(0,y,0X.07e+03 

a(0,0,z)l.07e+03 

yk 

P(0,0^^^ 


II 

i 


1* 


1 


a(x, 0,0) 

i 

3.99 

L 

a(0,y,0) 3.99 

a(0,0,z) 3.99 

P(0,0,v)' 

lllllllllllllllllllllllllii 

L ... 



F=0.47, G=3.41e+04, fi=0.0314. Ay/7=0.131 
THEORY: G(jB»l)=G. 73e+04 
THEORY: g (jE»l)=0. 0385 

a =499 at x =0.50, THEORY: a =420 at x =0.4G 
s s s s 

M^=0.1, c^(0,0)=0.4, c^(0,2)=0.14, c^(2,0)=0.12 


Figure 38. Results for N=130 of the third FEL. 
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V. CONCLUSION 


In this thesis, we have discussed the components and 
basic theory of a Free Electron Laser (EEL) in order to 
describe its behavior. 

We have explored, by means of simulations, three 
designs of Eree Electron Lasers in the amplifier 
configuration, using different parameters to establish a 
desired power output and optical wavelength. It has been 
shown that undulator taper has an important effect 
increasing the single pass extraction. 

Eor the low power (~100kW) EEL examples, extraction 
between 1% and 2% was found. Even after optimizing the taper 
parameters, the extraction did not improve much. We tried to 
obtain higher extraction by changing electron beam 
parameters, such as pulse duration, bunch charge, and 
emittance, but found no significant improvement. We conclude 
that these three parameters did not have a large effect on 
the final extraction of the system. 

In higher power (~3MW) FEL examples, it was found to be 
easier to obtain higher extraction. Eor high power EEL 
examples, we have obtained extractions above 3%, compared to 
1% or 2% in low power devices. 

Euture investigation should focus on taper parameters 
and undulator characteristics, as an important way to 
improve EEL designs. 
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